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Contrary to popular beliefs, it is possible to explain Baryonic asymmetry of 
the Universe inside the Standard Model, provided inflation ended into a bro- 
ken phase below the electroweak transition. Two important ingredients of the 
^ solution are multiquanta "Higgs bags", containing W,Z and top quarks, as well 

^ as sphaleron transitions happening inside these bags. Together, they provide 

I I baryon number violation at the level 10~^ -*. Our recent calculations show that 

CP violation (due to the usual CKM matrix of quark masses in the 4-th order) 
leads to top-antitop population difference in these bags of about 10~^. (The 
^ numbers mentioned are not yet optimized and simply follow a choice made by 

'"^ some numerical simulations of the bosonic fields we used as a reference point.) 
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1. Cold electroweak sceneirio 

It is a great honor to give a talk at this inaugural meeting of the Kobayashi- 
Maskawa Institute. Moreover, this invitation came at the moment whet I 
can report new exciting applications of the celebrated CKM matrix, to 
its originally intended purpose - explaining (perhaps) the most important 
CP-odd effect in our Universe, its baryonic asymmetry. 

The question how it was produced is among the most difficult open 
questions of physics and cosmology. The observed effect is usually expressed 
as the ratio of the baryon density to that of the photons UB/n^ ~ 10~^°. 
Sakharov"'^ had formulated three famous necessary conditions: the (i) baryon 
number and (ii) the CP violation, with (iii) obligatory deviations from the 
stationary ensembles such as the thermal equilibrium. Although all of them 
are present in the Standard Model (SM) and standard Big Bang cosmol- 
ogy, the baryon asymmetry which is produced by known CKM matrix is 
completely insufficient to solve this puzzle. 
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Significant efforts has been made to solve it using hypothetical "beyond 
the standard model" scenarios. In particular, possible large CP violating 
processes in the neutrino or in the supersymmetric sectors are possible. 
Although there are interesting proposals along those lines, they remain 
orders of magnitude away from feasible experimental tests. 

An alternative we will discuss is the modification of the standard cosmol- 
ogy. The standard Big Bang scenario predicts adiabatically slow crossing 
of the electroweak phase transition, leading to extremely small deviations 
from equilibrium. The so called "hybrid" or "cold" scenario^"^ solve this 
difficulty by combining the end of the inflation era with the establishment 
of the electroweak broken phase. Since there is no space here to discuss it in 
detail, let me simply enumerate the main points of the emerging scenario: 
only few new points will be elaborated below. 

• Large deviations from equilibrium, 0(1), with strongly oscillating 
Higgs and other fields. 

• Formation of relatively long-lived "hot spots" inside which the 
Higgs VEV is small.'^"^ 

• Topologically nontrivial fluctuations of the gauge field lead to 
baryon number violation^"^ at overall level of O{10~^) 

• These transitions take place only inside the "hot spots" .^"^ 

• The flrst nonscalar quanta produced are those with large mass (that 
is, stronger coupled to Higgs), namely the top quarks/antiquarks 
and the gauge bosons. 

• Tops and gauge bosons are collected into the "hot spots", which 
are then mechanically balanced. These objects are identified^^ with 
the (non-topological) solitons called the W-Z-top bags found inde- 
pendently^^ 

• Baryon number violation events are identified with Carter- 
Ostrovsky-Shuryak (COS) sphalerons with the size tuned to those 
of the hot spots 

• The B violation processes, described by the well known 12-fermion 
't Hooft operator, can occure as various subprocesses, n — )■ 12 — 
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n,n = 0..12 . The most probable was argued""^^ to be the "top recy- 
cling ttt — >■ 9 which is converting the energy of the tops or antitops 
already present in the bag into that of the gauge field, which helps 
passing high barrier separating different topologies 

• If so, the asymmetry of top and anti-top populations in the bag 
leads to different production rate of the baryon and antibaryon 
numbers^° 

• Clear "time arrow" of the process is given by the fact that a top 
quark, after a weak decay into lighter quarks, simply diffuse away 
from the bag since all quarks lighter than the top cannot be bound 

to the bag^° 

• The usual CKM mechanism of the CP violation naturally produces 
the top-antitop asymmetry via interferences of various outgoing 
paths of these light quarks, in the 4-th order in weak interaction^^ 



2. The multi-quark bags 

Being a scalar, the Higgs generates universal attraction between all kinds of 
particles. Furthermore, the strength of the attraction is proportional to their 
total mass, similar to the gravity interacting with the total energy. Gravity, 
feeble as it is, holds together planets, stars and even create black holes. 
Unlike vector forces induced by electric, weak or color charges, gravity and 
scalar exchanges are exempt from "screening" and thus their weak coupling 
can be compensated by a large number N of participating particles. And 
yet, unlike gravity, the Higgs boson is neither massless, nor even particularly 
light in comparison to W, Z or t. So, are there multi-quanta states based 
on the Higgs attraction? 

Another instructive analogy is provided by the nuclear physics. Think 
of a (much-simplified) Walecka model, in which the nuclear forces can 
be approximately described by the a, "the Higgs boson of the nuclear 
physics", and to meson exchanges. Because of similarity of masses ^ 
600 MeV, ~ 770 MeV, as well as couplings, the sigma-induced attrac- 
tion is nearly exactly canceled by the omega-induced repulsion. Their sum 
is an order of magnitude smaller than one would get from scalar and vector 
components taken separately. 

Can the situation at electroweak scale similar? Perhaps the Standard 
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Model is just a low energy effective Lagrangian, hiding some deeper physics 
behind its simplistic scalar Higgs. We are not aware of any particular model 
which suggests a vector companion to Higgs with a similarly small mass 
O{100GeV) mass. For example, the "techni-p" is predicted to be at the 
scale Mp ^ 1 TeV. Thus, unlike in the nuclear physics, one is not expecting 
the scalar-vector cancellation. 

The interest in the issue of "top bags" originated from the question 
whether a sufficiently heavy SM-typc fcrmion should actually exist as a 
bag state, depleting the Higgs VEV around itself. Although classically this 
seemed to be possible, it was shown in refs^^~^^ that quantum (one loop) ef- 
fects destabilize such bags, except at so large coupling at which the Yukawa 
theory itself becomes apparently sick, with an instability of its ground state. 
The issue rest dormant for some time till Nielsen and Froggatt^** suggested 
to look at the first magic number, 12 tops+antitops corresponding to the 
maximal occupancy of the lowest / = 0, j = 1/2 orbital, with 3 colors 
and 2 from t + i. Using simple formulae from atomic physics these authors 
suggested that such system forms a deeply-bound state. In ref.^^ we have 
checked this claim and found that, unfortunately, this is not the case. While 
for a massless Higgs there are indeed weakly bound states of 12 tops, they 
disappear way below the realistic Higgs mass. Further variational improve- 
ment of the binding conditions for the 12-quark system^^ confirmed that 
12 tops are unbound for Higgs mass mn > Mc(12) ~ bOGeV. 

Assuming spherical symmetry, the Higgs energy reads 



where v = 246 GeV ,m|j = 2Xv'^ is the Higgs mass, taken to be a round 
number''* rnn = 100 GeV in the original papers (wc also use; units of 
100 GeV throughout this paper). Consider now the addition of a conserved 
(during the time scale we are interested in) N particles (fermions or bosons), 
couple strongly to the Higgs field which could be strongly distorted. We 
adopt a mean-field approximation, in which all the particles are described 
by the same wave functions in the background of the Higgs field. Corrections 
to this mean-field description, such as, many-body, recoil and retardation 
of the Higgs field are expected to be suppressed by factors v/m, mH/m and 
1/N. In the semiclassical approximation, the total energy of the system will 
be given by E^i = E Higgs + X^a'^a^a' where {sa} is the spectrum of the 




(1) 



''If the ATLAS/CMS peaJc in diphoton will become a real Higgs mass, then it is about 
125 GeV. 
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corresponding field in the Higgs background. Ua is the occupation number 
of each state and N = Ua is the total, conserved, particle number. 

In the Higgs vacuum, i.e. (j){r) = 1, the state of lowest energy with N 
particles has total energy NM. However, in the background of a non-trivial 
Higgs field there are two competing effects. On the one hand, the gradient 
and potential terms increase the energy but, on the other hand, there might 
be some bound states levels with energy < < M which can allocate 
the quanta, lowering the energy of the system of particles at the expense 
of creating such distortion. 

Let us start by a crude estimate of the the order of magnitude of N for 
which such bags may exist. If we were to deplete a certain large volume of 
the Higgs VEV (surface/kinetic terms neglected for now), it would require 

2 

an energy Vbag^^v"^. For a bag of radius, say, R- 100GeV=4, this energy 
is about 20 TeV. Thus, if the lowest W-boson energy level has a binding 
energy of the order of 30GeT^ per or Z, an order of 0(1000) of them 
would be needed to compensate for the bag energy and obtain some binding. 
The top quarks are heavier and may get much larger binding, so one might 
naively think that less of them would suffice: but Pauli exclusion principle 
makes it more delicate. 

Consider the propagation of VF-bosons in an external Higgs field 



Let us study these equations of motion in the usual electric (e) , longitudinal 



(1) and magnetic (m) basis W^^'^'"^'' = Yj''^'"^^f(^g^i^^){r)/r, where Yj^ are 



spherical harmonic vectors and fe{r),fi{i') and fmir) are the radial wave 
functions for each mode. In a static, spherically symmetric, background 
the last term in (2) vanishes for the magnetic mode, leading to the simple 
Klein-Gordon equation: but in this case j > 1. For others one can start with 
smaller j = without a centrifugal potential, but the Laplacian mixes the 
electro-longitudinal modes, leading to the set of coupled equations (see the 
original paper) in which the last term in (2) becomes large and positive in 
the region where the Higgs field approaches zero, effectively repelling the 
longitudinal modes from the bag. (Note that massless gauge fields have no 
longitudinal degree of freedom at all.) As a consequence, even j = mode 
is pushed above that j = 1 magnetic modes, which is thus the lowest. In 
order to find bag solutions for finite N, we adopted a variational approach 





December 13, 2011 1:21 WSPC - Proceedings Trim Size: 9in x 6in shuryak proc KMI 



6 




2 4 6 8 10 12 14 16 18 20 2 4 6 8 10 

w R 



Fig. 1. (a) The energy of the lowest magnetic mode (j = 1) i?/(100 GeY for a Gaussian 
bag with a = 1.3 versus its width. The numerical results (points) are in good accordance 
with the analytical result, see.^* (b) Dependence of fermionic bound state levels energy 
£ in units of 100 GeV on the size of the bag, expressed as the parameter K in units 
1/100 GeV"^ for a bag with r\ = \. The color coding of the levels and their quantum 
numbers are listed in Table I. 



and took as a trial function for the Higgs, e.g. the Gaussian profile 

(\)[r) = 1 - a exp {-r^ jv?), (3) 

with two parameters, a and w describing its depth and the width, respec- 
tively. Solving the W-boson magnetic equation in this Higgs background is 
rather straightforward 

In Fig. 1(a) we show some results for a bag with a = 1.3. It is now 
relatively simple thing to vary the shape and reach the minimum of the 
energy of the system (still in the spherically symmetric ansatz.) 

Now we consider a system of N heavy fermions interacting with a back- 
ground Higgs field, the standard notations for Dirac spinors in spherical 
coordinates are 

where ^jim are spherical 2-component spinors and we take normalization 
/ dr {F^ + G"^) = 1. The so-called Dirac parameter k is defined as 

-(l + l) for j = 1 + 1/2 
forj = Z-l/2 



and runs over all nonzero integers, being positive for anti-parallel spin and 
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negative for parallel spin. The Dirac's equation reads 

{e-m(j))F = -G' + {K/r)G 
{e + m(j))G = F' + {n/r) F 



(6) 



The form of these equations presumes that the eigenvalue e is positive. A 
negative eigenvalue would correspond to a state in the lower fermion con- 
tinuum. If so, a charge conjugation transformation turns it into a positive 
eigenvalue for an antifermion. The Higgs equation of motion reads 

^2 {N - l)m F2 - G^ 



(7) 



Note that in the r — > oo limit, the source term in the right hand side, as 
well as the additional (p'/r term from the Laplacian can be neglected and 
the equation becomes the usual equation for a ID kink. 

The spectrum for Dirac's equation is foTind numerically and examples 
of the levels are shown in Fig.l. The Table shows magic numbers and the 
order in which levels are populated. Some levels are also shown in Fig. 1. 
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Attempting to find a minimum of the total energy, tops and Higgs, we 
have found that the ratio of the top-to-Hoggs masses are simply not large 
enough to stabilize bags by themselves. So tops can only exist inside the W- 
bags. Free (or weakly bound) top quarks are much heavier than W bosons 
and thus decay into another quark and the W. In the Higgs bag, however, 
we found that two lowest top levels are below the lowest gauge boson ones; 
so up to 36 t + f in such bags may live much longer, given by next order 
decays into three fermions, like in the usual beta decays. 

3. The sphalerons and rates 

In the broken phase the electroweak sphalerons had been found by Manton 
et al: their mass is about 14 TeV and thus the tunneling rate is prohibitively 

small. In pure gauge sector finding the sphaleron solution has been pre- 
cluded by the fact that classical gauge theory has no scale and thus energy 
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has no minimum. It has been surpassed in^^ by requiring minimization un- 
der two conditions: fixed Chern-Symons number and mean-square radius 
P- 

The profile of the the magnetic field in COS configuration is given by 
the following spherically symmetric expression 

^4 



B\r) = ^.p^ (8) 

This form does indeed fit well the numerically found shapes of the at the 
"sphaleron moment" , a maximum of magnetic field. Fitted radius, yielding 
mp « 3.9 which corresponds to the total energy of the COS sphaleron 

Etot = ^n^/g^P « 2 TeV (9) 

which is 7 times less than the KM sphaleron mass: it makes a huge diff'erence 

for he rate. 

Combining all the factors we find that numerical value of preexponent 
and exponent nearly cancels out, leaving crudely 

r/Fm^ ~ 10"^ (10) 

with accuracy say an order of magnitude or so. With that accuracy it 
agrees with the results of the simulations which also finds that the number 
of sphaleron transitions per spot is indeed about several percents. 

What should happen after the sphaleron moment? Sphaleron decay is a 
classical downhill rolling of the classical (high amplitude) gauge field, from 
the (sphaleron) top into the next classical vacuum . This process was ex- 
tensively studied numerically for the broken-phase sphaleron. Remarkably, 
an analytic solution of the time-dependent explosion of COS sphaleron has 
also been found in the COS paper. The late time profile of the energy 
density of the expanding "empty" shell. 

Stt / 1 \ 

'"^(^'^)^g^( l + (r-t)VpO ^''^ 
Comparing the explosion of COS sphaleron with numerical data one can 
see both the similarities and the differences between them. Indeed, there is 
an empty shell formation at some time. However the inside of the shell does 
not remain empty: in fact the topology and magnetic field have a secondary 
peak (of smaller magnitude). Qualitatively it is easy to see why it happens. 
The COS sphaleron is a solution exploding in zero Higgs background, with 
massless gauge fields at infinity. In the numerical simulations we discuss 
such explosion happens inside the finite-size cavity. As the gauge bosons 
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of the expanding shell hit the walls of the no-Higgs spot, they are massive 
outside. With some probability they 

Since the original numerical simulations have included the gauge fields 
but ignorcicl ferniions, wc; have to discuss first, at quite qualitative level, 
what their effects can be. relative to that of the gauge fields. Those tops 
would be added to the metastable bubbles of the symmetric phase, the 
no-Higgs spots, like the W discussed above. 

The well known Adler-Bell-Jackiw anomaly require that a change in 
gauge field topology by AQ ± 1 must be accompanied by a corresponding 
change in baryon and lepton numbers, B and L. More specifically, such 
topologically noritrivial fluctuation can thus be viewed as a "t'Hooft opera- 
tor" with 12 fermionic legs. Particular fermions depend on orientation of the 
gauge fields in the electroweak SU(2): since we are interested in utilization 
of top quarks, we will assume it to be "up" . In such case the produced set 
contains trthtgCrCbCgUrUhUgT^i, e , where r,b,g are quark colors, to which 
we refer below as the — > 12 reaction. Of course, in matter with a nonzero 
fermion density many more reactions of the type n — > (12 — n) are allowed, 
with n (anti)fermions captured from the initial state. 

The classical solution describing the expansion stage at i > has been 
worked out for COS sphaleron explosion, and for the "compression stage" 
at f < one can use the same solution with a time reversed. At very early 
time or very late times t — >■ ±oo the classical field become weak and describe 
convergent /divergent spherical waves, which are nothing but certain num- 
ber of colliding gauge bosons. Fermions of the theory should also be treated 
accordingly. Large semiclassical parameter - sphaleron energy over temper- 
ature - parametrically leads to the assumption that total bosonic energy 
is much larger than that of the fermions, so one usually ignores backreac- 
tion and consider Dirac eqn for fermions in a given gauge background. For 
KM sphaleron and effective T we discuss, this parameter would be ~ 70, 
which is indeed large compared to 12 fermions. However in the case of COS 
sphaleron we are going to use the number is about ~ 10, comparable to the 
number of fermions produced. It implies that backreaction from fermions 
to bosons is very important. 

The only (analytic) solution to Dirac eqn of the "expansion stage" was 
obtained in,^'' it describes motion from the COS sphaleron zero mode (at 
t = 0) all the way to large t — )■ +00 physical outgoing fermions, with an- 
alytically calculated momentum distribution. A new clement pointed out 
in^*" is that its time-reflection can also describe the "compression stage", 
in which free fermions with the negative energy are captured by a conver- 
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gent spherical wave of gauge field at t — )• — oo, ending at the zero energy 
sphaleron zero mode at t = 0. 

This implies that energy of the initial fermions can be incorporated and 
used in the sphaleron transition. The optimum way to generate sphaleron 
transition tm'ns out to be 3 initial top quarks'' considering the 3 — > 9 
fermion process instead of the original — >■ 12 one. The 3 — )■ 9 fermion 
process saves a lot of energy, as in it the initial top quark energy can be 
completely transferred from the "sphaleron doorway state" to the gauge 
field. Estimates show that it increases the sphaleron rate by about one 
order of magnitude, compared to pure gauge calculation. 



4. The CP cisymmetry 

The first attempts to estimate magnitude of CP violation in cold elec- 

troweak cosmology has been made by Smit, Tranberg and collaborators. ^'^""^ 
Their strategy has been to derive some local effective CP-odd Lagrangian by 
integrating out quarks, and than include this Lagrangian in their real-time 
bosonic numerical simulations. The estimated magnitude of the CP-odd 
effects ranges from ns/ny ^ 10~^..10~^'',^ which reignites hopes that this 
scenario can provide the observed magnitude of the baryon asymmetry in 
Universe. However, there are many unanswered questions about the ac- 
curacy of these estimates. One of them-^'' is that the effective Lagrangian 
derived with specific scale of the loop momenta, e.g. p ^ rric, can only be 
used for field configurations at a scale softer than this loop scale: and the 
"hot spots" in numerical simulations obviously do not fit this condition. But 
in practice even more important is the following unanswered generic ques- 
tion: why should a very complicated operator (containing 4-epsilon symbol 
convoluted with 4 gauge field potentials and one field strength) averaged 
over very complicated field configurations (obtained only numerically) have 
nonzero average at all? We were thinking aboiit some model fields (e.g. 
sphalerons) in which these operators have nonzero values, but were not 
able to find any convincing examples. Since the calculation is numerical, it 
would be desirable to have some parametric estimate of the effect, in par- 
ticular to know what sign the effect should have and at least some bound 
on it from below. These goals are reached in the last paper. 

For a top quark starting at the position ri , the escaping amplitude has 



"In order to satisfy Pauli principle and fit into the same sphaleron zero mode, colors of 
the 3 quarks of each flavor should all be different. 
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the form. 

At{ri) = jYW-V^PtS^{rur2h^W+VS^{r2,r,)d^r2, (12) 
whereas the anti-top has a C-reflected expression 

Atiri) = Jj-W+V^PtS!^in,r2h^W-VS^{r2,r,)dW (13) 

Here V is the CKM matrix, S the quark propagators, their index u, d etc 
denotes the up and/or down quark flavors and Pt denote the flavor matrix 
making the initial projection on the top quark. 

The probabhhty of a top quark escaping from n is then given by the 
integral over all positions and sum over all intermediate and flnal states / 
of the squared amplitude 

Probt{n) = J^AlAt=Tr J ^^,^4^*= J d^r,d''r2d^r3Tr [Ptj''W-V^ S^{n,r2) 

rW+VSt{r2,rc)S^\ro,r3) rW-V^ S^\rs,nh^W+V . 

Note that the interference terms between different paths are of the 4-th 
order in weak interaction, and thus have 4 CKM matrices, as indeed is 
needed for the CP violation effects. Four positions of the points at which 
the interactions take place, as well as particular quark flavor in the interme- 
diate line, are summed over. Writing the amplitude squared of the process, 
one includes the unitarity cut (the vertical line in Fig. 2) to the right of 
which one, as usual, flnds the conjugated image of the process in opposite 
direction. 

In between these four points the flavor of the quark remains unchanged. 
Quark wave functions (we keep in mind Z = or s-wave ones only, thus 
points are only indicated by their radial distance from the bag) are differ- 
ent for each flavor, because of different Yukawa couplings to Higgs profile. 
Semiclassically the phase is approximated 

Si2 = exp[i J pix)dx] « exp[i j {E - -^)dx], (14) 

where E is the quark energy, and the approximation implies that all lower 
quark flavors are light in respect to r / E, so that the flavor-dependent phase 
(stemming from the second term in the bracket) is still smaller than 1 

Let us follow the flavor part of the amplitude, which distinguishes be- 
tween quarks and anti-quarks. The 4-th order process we outlined in the 
preceeding section corresponds to the trace of the following matrix product 

Mt = Tr{Pt * V * Sf^"-'" * V+ * S^i * V * S-^^" * ^^). (15) 



December 13, 2011 1:21 WSPC - Proceedings Trim Size: 9in x 6in shuryak proc KMI 



12 





b,S t 



W 



Fig. 2. (Color online) Schematic shape of the 4 order process involving only quarks 
of the 2nd and 3ed generations. The shaded objects on the left and right represent the 
Higgs bag with strong gauge fields (indicated by W in the figure) inside. The vertical line 
is the unitarity cut. Four black dots indicate 4 points, indicated with ri,i = 1..4 where 
W quanta are interacting with the quark, changing it from upper (t=top) to lower (fe, s) 
components. 



where S are quark propagators, the lower indices specify their initial and 
final points, the upper subscript remind us that those are for up or down 
quark components. Pt is the projector requiring that we start (and end the 
loop) in the bag, with a top quark. We also define the amplitude for the 
top-antiquarks 

Mi = Tr{Pt * V* * Sf^"-'" * V'^ * S^^ * V* * Sl^" * V^) (16) 

which we subtract from Mt, as the effect we evaluate is the difference in 
the top-antitop population inside the bag. The difference gets CP-odd as 
seen from its dependence on the CP-odd phase S 

Mt — Ml = 2iJ (5*23 — *S'23) (— S"!]^ * — S'i^i * 5*^2 + * 

+3^2 * — * Si^i + SI2 * S^i) (17) 

J = sm{S) sin(6'i2) sin(6'i3) sin(6'23 cos(6l23) cos(6'i2) cos^(6'i3) cos(6'23) 

The remaining combination of propagators, organized in two brackets, 
needs to be studied fiurther. Note first that the propagators in the range 
2-3 (through the unitarity cut) factor out and that one may ignore the top 
quarks there. Note further, that if the u, c quarks would have the same mass, 
the first bracket would vanish: this is in agreement with general arguments 
that any degenerate quarks should always nullify the CP-odd effects, as the 
CP odd phase can be rotated away already in the CKM matrix itself. 

The last bracket in (17) contains interferences of different down quark 
species: note that there are 6 terms, 3 with plus and 3 with minus. Each 
propagator, as already noticed in the preceeding section, has only small 
corrections coming from the quark masses. Large terms which are flavor- 
independent always cancel out, in both brackets in the expression above. 
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Let us look at only the terms which contain the heaviest b quark in the 
last bracket, using the propagators in the form 5*^'^ ^ exp ^zti^rij^ where 
± refers to different signs in the amplitude and conjugated amplitude and 
rtj = Tj — rt. Note that the sign of the phase between points r2 and rs can 
be positive or negative as it results from a subtraction of the positive phase 
from rs to the cut Vc with the negative phase form the cut Vc to r2- Terms 
containing odd powers in r23 therefore should vanish in the integral: and 
the lowest term we have is quadratic. Considering all phases to be small 
due to 1/E and using the mass hierarchy nii, ^ ttIs ^ we pick up the 
leading contribution of the last bracket in (17) which has and the 4-th 
power in the last bracket, the 6-th order in the phase shift in total: 

g4^5 V ; 

Note that all distances in this expression are defined to be positive and 
the sign in the last bracket is plus, so unlike all the previous orders in the 
phase expansion, at this order we have sign-definite answer with no more 
cancellations possible. This point is the central one in this work. We further 
see that this expression grows for large r's, which are to be integrated over. 
Of course as we expanded the exponent in the phases, they have to be such 
that these phases are smaller than 1. 

Let us start with a "naive" estimate, which assumes that E in the 
formulae is given by the top quark mass E rrit = 173 GeV. As for the 
field strength, naively one may take all four interaction points inside the 
bags, where the amplitude of the W is the strongest. If so, all distances r^j 
are of the order of the bag size Rbag ~ l/mw However, if this is the case, 
all the phases are so small that the resulting CP asymmetry is about 10 
orders of magnitude smaller than needed. However, the initial top quarks 
are bound in the bag, so light quarks can be propagating at the energy 
much smaller than the top mass. The smallest possible scale is fixed by the 
weak interaction of quarks with the electroweak plasma outside the bag, 
known as the screening mass ^ g-wT, which is few GeV. This is the natural 
scale to take: thus we will from now consider E ~ mt, in the following. 
Another improvement one may try is to consider locations of some points 
outside the bag, selecting rij as large as possible. 

Probt_,{n) ^ J^^I^IP^ J dr2drs2r2rlf(nrf(r2).f{r3). (19) 

Considering a radial bag of Nw weak bosons having an exponential profile 
with the usual W mass in the broken phase: W{r) = ^^^"^ Q-m„r^ 
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get that the probability of a top-minus-antitop quark escaping 
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Sep = tJN^ 
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10 



(20) 
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In the latter formula we made use of the lifetime of the bag denoted r/m^y , 
with r 6 to bound the time integral over Xi. 

The main lesson we got from this study is that the scales of both the 
quark energy E and their traveling distances in the loop amplitudes 
should be tuned individually, to maximize the effect. The main limitation 
come from the conditions of quark rescattering in the plasma (the screan- 
ing masses) and the conditions that all phases Si should not be large, as 
well as the limitations coming from the W field strength and correlation 
length. Another lesson is that in order to prevent cancellations between 
different flavors, one has to expand all the results till sign-definite answer 
is guaranteed. 

The probability to find 3 antitops is actually proportionnal to (1 + 
Scp)^ w 1 + 3Scp, while it is (1 — 3Scp) for tops: it gives factor 3. Another 
factor 3 appears because of the fact that each sphalcron event creates 3 
units of baryon number, not one. Together with baryon asymmetry (time 
integrated) sphaleron rates of 10~^ and 3 * 3 * ^c-p we arrive to our final 
estimate AB ^ 10~^^^^ where one order of magnitude stands for our (per- 
haps optimistic) errors due to mimcrical factors ignored in the estimates. 
We conclude that it is clearly in the same ballpark as the observed baryonic 
asymmetry of the Universe. Clearly, numerical factors can be detailed later, 
and the parameters of the cosmological model can be better tuned to get 
closer the right value. 

Last but not least is the issue of the sign of the asymmetry. Our formula 
(18) has definite (positive) sign, that is to say more top quark escape the 
bag (note that the time direction is important, quarks arc first created in 
the bag, then have more probability to escape). More antitops remain in 
the bags, with more likely to be "recycled" by the sphalerons: this produces 
more baryons than anti-baryons. Apparently we got the right sign for the 
baryon asymmetry. 
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